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DISTURBANCE OF OPERATION IN QUANTUM 
ESTIMATION FOR THE GAUSSIAN P-FUNCTION 

Yoshiyuki Tsuda* , Keiji Matsumoto* and Masahito Hayashi** 

For the quantum Gaussian state family, Hayashi proposed a quantum mechan- 
ical operation using beam splitters to estimate the location and scale parameters of 
the P-function, and he showed that it is asymptotically optimal. In this paper, we 
^>f-^ ' analyze the effect of disturbance of his operation caused by the randomness of the 

I transparency of the beam splitters. It is shown that even if the variance of the ran- 

. dom transparency is small, Hayashi's estimators are improper in a sense that they 

CNj ' are biased and asymptotically inconsistent. In such a case, we propose to stop the 

5_j operation and correct the biases of estimators. 

• Key words and phrases: coherent state, counting measurement, quantum Gaussian 

state, heterodyne measurement. 

00 ■ 



(N 



1. Introduction 

^ . The quantum estimation is an application of the statistical estimation theory to 

the quantum mechanics. The general theory on the quantum mechanics is described by 
the linear algebra on a Hilbert space, and it might require some specific setup in order 
' to consider the quantum estimation. However, our present problem on the Gaussian 

. P-function is a special case in which only the basic analytical methods are used. This 

' paper is unfolding our problem without linear algebraic premises. 

O ' m 

' The quantum estimation theory was first considered by Helstrom [2]. The main 

, problem is to find a optimal scheme to estimate unknown parameter of the quantum state. 

II ' Here, the scheme of estimation is composed of (i) physical operation and measurement 

, in order to obtain data and (ii) the estimator based on the obtained data. Let p denote 

' a state of the system to be observed. Let M denote a measurement on the system. If 

^ , one carries out a measurement M, then he will obtain data x as an observation of a 

O^' random variable X whose probability distribution is determined by p and M . Suppose 

that we know the true state is an element of the set S := {pg} parameterized hy 9 € Q, 
and also suppose that we can select a measurement in the set M = {Mx} parameterized 
by A e A. Based on the obtained data x, one will estimate S by an estimator 9. 
' In the standard statistical estimation problem, a statistician is allowed to select only 

the estimator 9. One features of the quantum estimation is that one could improve 
the accuracy of the estimation by interaction between samples. In the present paper, 
we consider the effect of disturbance of physical operation in a problem of quantum 
estimation for the Gaussian P-function model. 

The models of Gaussian P-function are typical and important examples of the quan- 
tum estimation and it can be applied to the optical experiments. Our model has two- 
dimensional location parameters and one scale parameter. For the estimation problem 
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of the location parameters of the Gaussian P-function, Yuen and Lax and Holevo 
^ have given a lower bound of the variance of unbiased estimators and an optimal 
estimation scheme that attains the bound. For the problem of the location and scale pa- 
rameters, Hayashi ^ has given a lower bound of the variance of unbiased estimators and 
an optimal estimation scheme that asymptotically attains the bound. His result gave the 
first practical example of the advantage of interaction of samples to the quantum esti- 
mation of the Gaussian model concerning the statistical estimation theory. We consider 
the model of the location and scale parameters and analyze the effect of disturbance of 
the asymptotically optimal scheme given by Hayashi . He pointed out that his scheme 
is realized in the optical experiment by the physical operation using beam splitters. For 
example, if we have given n = 2™ samples, then an experimenter uses beam splitters 
whose transparency is 1/2. However, in practical, the transparency of prepared beam 
splitters is not always exactly the ideal quantity (1/2), and it is often slightly different 
from the ideal one. If the transparency is not correct, then the physical operation pro- 
posed by Hayashi IT] is disturbed and it might affect the estimation. We analyze the 
effect of the disturbance to the estimation, and it turns out that the estimator is not 
even asymptotically consistent to the true parameter. In such a case, a naive scheme not 
using interactions between samples is better than that of Hayashi ^ since the distur- 
bance could be avoided. In order to avoid the inconsistency of Hayashi pP's estimator, 
we propose a new scheme to estimate location and scale parameters by correcting that 
of Hayashi [T|. We will see that our scheme gives an asymptotically consistent estimator, 
and it is better than the naive estimator. 



2. Problem 

When we discuss a general problem of quantum estimation, we need to consider 
any measurement M which may be carried out to the quantum system. However, our 
interest in this paper is devoted to the problem of disturbance in the experiment proposed 
by Hayashi jlj. His optimal method uses only two kinds of measurements called the 
'heterodyne measurement' and the 'counting measurement' and the state is denoted 
by a probabihstic superposition of the coherent states. It enables us to describe the 
problem by a Bayesian model in which the prior distribution of parameters of basic states 
means the state, and selecting some of random variables corresponds to the selecting 
measurements. This is why the linear algebraic preparation is not needed. 

Let {Aj, Bj) denote the complex amplitude Aj +iBj {i — ^J—l) of a coherent state, 
{Xj^Yj) denote the observed amplitude Xj + lYj by the heterodyne measurement and 
Zj denote the counted numbers by the counting measurement. Suppose that Wj := 
{Aj, Bj, Xj, Yj, Zj) (j = 1, n) are random variables. Suppose that we can always not 
observe Aj or Bj for any j, and we can selectively observe {Xj, Yj) or Zj for each j, so 
we do not observe Zj if we select to observe {Xj , Yj ) and we do not observe {Xj ,Yj) if 
we select Zj. Let S C {1, n} be a set of indices of which we observe {Xj, Yj), hence 
we observe if fc ^ S. Let {xj,yj) {j G S) and Zk {k ^ S) be observations of {Xj,Yj) 
and Zj. For each j g S, we assume that Xj and Yj are continuous random variables 
on the set R of real numbers, and, for each k ^ S Zk is a discrete random variable on 
the set No of non-negative integers. We also assume that the joint probability (density) 
function of {Xj, Yj) {j E S) and Zk {k ^ S) is denoted as 



E 



j-j / eM~{^,~A,)'-{y,-B,)') \ j-j / cM-Al-Bl){Al + Bl)^ 

j£S V / k^S ^ 



where the expectation is taken for Aj and Bj {j ~ 1, ...,n). 



REMARK 1. Suppose that Ai, An, Bi, Bn are independently distributed. 
For j e S, Aj and Bj are distributed according to the normal distributions N[0,v/2) 
and N{'ii, v 12) of means 9 and rj and variance i//2, respectively. For k ^ S, Ak and are 
distributed according to A^(0, J//2). Then, all selected variables {Xj,Yj) {j g S) and 
{d ^ S) are independently distributed. The distributions of and Yj are N{9, (i/+l)/2) 
and N{r], {i/ + l)/2), respectively, and that of is the geometric distribution Geo{v) of 
mean that has the probability function 



Not only we can select observed variables, but we can also select pairs of variables 
{Aj,B.j) and {Ak,Bk) and transform them as 



where r e i? is arbitrarily selected. We call this transformation g^' ■ Note that the 
transformations has the group structure in the algebraic sense, and we write gii^ ''^^^ gii^ ''^^'^ 
to denote that we first carry out gii'''^^^ and then next do g^^'^'^^ ■ 

REMARK 2. Suppose that Aj, Bj,Ak, are independently and identically dis- 
tributed. Let Aj and Ak are distributed according to N{9,v/2) and Bj and Bk are 
distributed according to N{ri,v/2). If we carry out the transformation gT^'^\ then the 
distributions of Aj, Bj, Ak, Bk are given by 



respectively. 

Suppose that A^ , A" , Bi, Bn are independently distributed, and Aj (j — 
l,...,n) is distributed according to N{9,v/2) and Bj {j = l,...,n) is distributed ac- 
cording to N{r],p/2), where 9, rj Cz R and v > are unknown. Under the rules of 
selecting observed variables and transformations, we consider to estimate 9, rj and v 
based on observed variables. 

3. Naive estimation Vifithout transformation 

In this section, we consider one of the most naive method to estimate unknown 
parameters, that is, we do not carry out any transformations and we let S := {1, n}. 
Then, from Remark Q Xi, Xn,Yi, ...,Yn are independently distributed, and Xj is 
distributed according to N{9, (i/+l)/2) and Yj is distributed according to N{ri, {i'+l)/2) 
for j = 1, n. Let 





iV(6'(cosT + sinT),j//2), N{r]{cosT + sinr), i^/2), 
N{9{- sinr + cost), iy/2), iV(77(- sinr + cost), iy/2). 
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be estimators for 9^ rj and v. Then, we have 

E{e) = 0, E{n) = V, E{D) = u, 
so that they are unbiased. Moreover, since the covariance matrix of {9, fj, v) is 
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we can see that they are asymptoticahy consistent. 

4. Hayashi's estimation using transformation 

In this section, we consider Hayashi ilj's method to estimate 0, t] and v. Note 

first that, by some transformations of Ai, An, Bi, Bn by G := g^n 9ti'''^\ 
we obtain random variables A[, A'^, B[, B'^. which are mutually independent and 
Ai,Bi, Aj, Bj (j > 2) are distributed according to 

(4.1) NiV^9,iy/2), N{^ij,v/2), 7V(o,i./2), iV(0,j//2),, 

respectively. 

REMARKS. For t = 1, ...,n- 1, let 

Tt := tan"4"^/2. 

Then, a transformation defined by 

Gi 

generates A[, A[^, B[, B'^^ which are independently distributed according to H4.1|l . 



If we draw the transformation gr''^'^ as 



iA,,B,) 

{Ak,Bk) 
then, Gi can be expressed as 



{Aj cos T + Ak sin t, Bj cos t + Bk sin r) 



-e- 



(4.2) 



iAi,Bi) 

(^2,^2) 
(^3,^3) 

{Ai,Bi) 



'Aj sin T + Ak cos r, — i?j sin r + _Bfe cos r) 



{A[,B[) 
{A'„B',) 
{A'„B',) 
{A',,B',) 



-e- 



7r/4 



-e- 



tai-i 2-1/2 



-e- 



tan-i 3-1/2 



for n — 4. 

The transformation generating A'l, A^, B'l, .... i?^ which are which are distributed 
according to H4.1|) is not unique as the follows. 

REMARK 4. Suppose that n = 2™ for a natural number m. For t = 1, to, let 

Jt '■= {(jj ky\ be the set of pairs of natural numbers j and k satisfying 



j = l (mod 2*), fc-j 



Then, a transformation defined by 



G2 



n 



■ n 



n U,k) 
y7r/4 



also generates A[, A[^, B[, B'^ which are independently distributed according to 

6U. 



In the same way as H4.2|l . G2 can be expressed as 

(^2,52) 
(^3,^3) 
(^4,^4) 



{A'„B',) 
{K,B',) 
for n = 2^ = 4. 

After such transformations obtaining A[, A'^^, B[, B'^, we let 5* :— {1} to select 
observed variables {Xi ,Yi), Z2, Zn. Let 
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t < 

7r/4 
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Xij^/n, f] := Yi/vn, 



1 " 



i=2 



then from Remark ^ we have 

so that they are unbiased. Moreover, since the covariance matrix of (0, ?}, is 

1^2 := 



1^+1 

2n 










2n 










n— 1 



we can see that they are asymptotically consistent and this Hayashi's estimators ^, 77, v 
dominates the naive estimators ^, 77, j/ in a sense that the difference V\ — V2 of the co- 
variance matrices is non-negative definite. 

Hayashi ^ also proved that his estimators are asymptotically optimal in the quan- 
tum mechanical setup which is more general than that of here. 



5. Noisy transformation 

The transformations G of Ai, An, Bi, Bn to A[, A'^, B[, B'^^ is closely 

(i k) 

related to physical operation of in quantum optics. Actually, each g)- corresponds 
to the interference of two light beams using a beam splitter of transparency cos^ r. For 
example, if an experimenter tries to realize the transformation G2, then he has to prepare 

(i k) 

beam sphtters of transparency 1/2, since G2 is constructed by g]^j^"s. 

However, in practical cases, it is difficult to prepare beam splitters of exactly the 
same transparency as the ideal quantity (1/2), and in many cases, each quantity is slightly 
different. Hence, we consider, when n = 2™ and we try to carry out transformation G2, 
but each r is independently and identically distributed according to A^(7r/4, elog2), and 
we calculate expectations and variances of Hayashi's estimators 6, fj, v in order to see the 
effect of the randomness of transformations. 



Let A' :— {A[, ...,A'^, B[, B'^). Then, note that Remark ^ shows that 
E{e) = ^E{E{X, I A)) = 
E{e^) = -^E{E{X2 I Aj) = E{B[), 

n 
n 



= -Eiyix^ I A)) 

n 

= -S(Af + 1/2) - EiQf, 
n 

v{f,) = ii?(y(ri I A)) 

n 

= ii?(i3f + 1/2) - 

n 



= 7;r^ E (^((4 + B',f + (4 + B'^f) - EiA- + Bff) 
^ ' i=2 

E {mA',+B'^nA',+B',f)) 

^ ' 2<j<k<n 

E {E{{A'^+B'^r)E{iAi+B',r)), 

^ ' 2<j<k<n 

hold, so that these quantities depend only on bivariate, at most, fourth moments of A, 
and we can obtain the following results. 

THEOREM 1. For any e>0, we have 

(5.1) E{e) = n-'^^e, E{fi) = n-'/^ij, 

(5.2) Vie) 



1 + z/ if' rP' n^^ — n^'- , 

(5.4) Ei,) ^ . + -l + 2^-^-2^-^n,-2^^n-,^^n^-^ 

^ ' ^ ' 2'^(-2 + 2^)(-l + n) ^ ' ' 

If e is positive and small enough, then we have 

(5.5) V{v) = 0{n-^') 
as n goes to oo. 



Af and B"-^' after the transformation Ilo-.fcjGJt 5^'/^ ' ' ■T\(]M)eJ^ y^/A 
way of G2- Then, for any natural numbers ta, tb, sa, sb, t and j — 2*^^ + 1, the 
recurrent equations 



PROOF. First, for j = 1, ...,n (= 2™) and < = 0, 1, ...,to, let Af and B^ denote 



,(t) 



which is on the 



yl^*~^^sinT) 



(t-i) , 



Sf"^) sinr)'^«), 



(t-l)rAQ{t-l)rB\ 



E{A'I 



E{A 

(t-l)rA -g(t-l)rB ^(t~l)sA Jj(t-1)SB 



— E{A^^^^^^'^ B^^^^^^^) E{A^*~^^^'^ B^*~^'^^^) 



hold. Hence, we have 1)5. l|l to (|5.3|l . Next, for any ta, rs, t, j = 2* ^ + 1, and for any 
natural number k satisfying 



(mod 2*), 



the recurrent equations 



£;((-A^*"^'sinT 



A*!*"^) cost) 



B 



sin T + B 



(t-i) 



cos ry^), 



E{A 



{t)rA r>{t)rB\ 



- E{A^,>''^B)^ 



hold. Hence, we have H5.4() . Finally, for any r^, r^, sa, ss, A = 1, m— j = 2* ^- 
and k — 2*+^^"^ + 1, the recurrent equations 



WrA Jj{t)rB ^(t)sA j^{t)sB 



E{A'^''^B 



4 



= E{{A';''' COST + sinT)'--^ (Bf cost 



x{-A['^^hmT + A'-''"' cosTy^{~B[ 



(t-i) 



sinr + Bf costY^), 



E{{A['+^-'^ cost + siuT)'-- 

sIut)'"^ 

x{-A['^^'^ siuT + Af^^^ cos tY^{-B[''^^ sinT + cost)'*^), 



k 

i 

= SinT + 4*+''-') cost)--- 



X (-sf+^-i) SinT + cos Tf- 



x(-Af siuT + cosT)''^(-Bi'"'^ SinT + B)"^'' cost)''^) 



□ 



hold. Hence, we have (|5.5|l . 

By this theorem, we can see that Hayashi's estimators 9, fj and i> are biased and 
asymptotically inconsistent. 



Let us compare v and v. Since v is biased, it is better to compare them by mean 
square errors, that is, 



M E{(p-vf) = V{v), 

If n = 2, then we have 
M = + 

M = ^^ + ;/+(2;/+l)(l-2-^^)(0^ + r;^) + ^"^^ ^"^ \ o^ + V^)- 

We can see that, if e or 6*^ + ry^ is small enough, then M < M. In case of 77 = 0, (e, 9) 
satisfying M — M are numerically calculated for v = 0.1, 1, 10, 100 (solid lines) and for 
u = 0,00 (dashed lines) and plotted in Figure ^ If the parameters are in the lower-left 
side of the lines, then M < M, so that Hayashi's estimator 1/ for v has advantage in the 
sense of mean square error. 



Figure [I 



6. Correcting Hayashi's method by stopping the transformation 

We have seen that Hayashi's estimators 9, ij and i> are not so good asymptotically 
if the transformation is noisy. Hence, in the case of n = 2™, we consider to stop the 
transformation G2 of Remark 01 at mo-th step. Let Wj := [A'j, B'j, Xj,Yj , Zj) {j = 
1, n) be random variables obtained by the transformation 

n n 9i''"^--- n 9^^-'^ 

of iA,,B,,X,,Y,,Z,) [j = 1, n). Let 



s 



:= {2'"«j + l I j==0,...,2™-'""-l} 



be a set of indices of which we observe {Xj,Yj), hence we observe li k ^ S. Let 



1 2™^'"o — 1 2 ^ 

1 I + ('2m _ n2m-mQ ^^^'''i ^ ^ 

k^s ^ ' jes 

a ^ , , 2™-2™« 



where 



_|_22)no+e(2+nio) _ 2"'+2™0+f (2+mo) _ g™0^ 

/3 := (-2 + 2')(-l + 2"')(-2" + 2"«). 
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Figure 1. Pairs of (e, 6) satisfying M = M for v = 0.1, 1, 10, 100 (solid lines) and for = 0, oo (dashed lines) in case of 



Then, by the same argument of the proof of Theorem ^ we have 



{2" - 1)2(2"'« - 2"=™") 



V{flmo) 
and, if mg is fixed, then 

as n goes to oo. 



2l+m — emo 
2l+m — erriQ 



2e+m(2<= - 2) 

{2" - 1)2(2'"" - 2'=™") 

2e+m(2e - 2) 



7. Numerical comparison and conclusion 

For n = 2^28,210,212 and for e = 0.0001, 9. 
for 6 = 0,2, 4, 6, 8, 10, ry = and v = 0,2, 4, 6, 8, 10, the sum of mean square errors 



, "mo: Vmo ^-^d i>rno coustructed, and. 



E^mo - of + iVmo - Vf + iymo - v)) 



are numerically calculated and their relative errors 

l-M™jMo 

are shown in Table 1. The Hayashi's estimators or the corrected estimators are better 
than the naive estimators if the relative error is positive. The corrected estimators 
^m/2j f\mli and I'm/2 are better than the Hayashi's estimators if the relative error of the 
corrected ones is larger. 

When the location parameter d is close to the origin or the scale parameter v is large, 
the effect of the randomness of the transformation is small and Hayashi's estimators has 
the good performance. However, when the location is far from the origin or the scale is 
small, the randomness of the transformations significantly influences Hayashi's estimators 
and its performance is inferior to the naive estimator. Our corrected estimators also loses 
the accuracy when the location is far from the origin or the scale is small, but such bad 
effects are relatively smaller than that of Hayashi's estimators. Indeed, in the range of 
our simulation, the corrected estimators are always better than the naive estimators (see 
Figure 12 for n = 2^^). 



Table 1 



Figure 12 
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Figure 2. The goodness of estimators for n = 2 



